Percolation model for nodal domains of chaotic wave functions 
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Nodal domains are regions where a function has definite 
sign. In it is conjectured that the distribution of nodal do- 
mains for quantum eigenfunctions of chaotic systems is uni- 
versal. We propose a percolation-like model for description of 
these nodal domains which permits to calculate all interest- 
ing quantities analytically, agrees well with numerical simula- 
tions, and due to the relation to percolation theory opens the 
way of deeper understanding of the structure of chaotic wave 
functions. 
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In a recent paper Smilansky et al. consider the 
following problem. Let '^{x,y) be a real eigenfunctions 
of a 2-dimensional quantum problem. The equation 
^'(x,?/) = determines a set of nodal lines which sep- 
arate nodal domains where '^{x,y) is of opposite signs. 
In 1^ it is argued that the distribution of the number of 
these regions for high excited states is (i) universal for 
integrable as well as for chaotic models but (ii) it clearly 
distinguishes between these two types of models. 

For chaotic (billiard) systems it is conjectured in ||] 
that this distribution coincides with the distribution of 
nodal domains for Gaussian random functions which are 
known to give a good description of wave functions of 
chaotic systems pi: 



C„vE'(^)(x,y), 



(1) 



where ^>m\x,y) — J|m| (fcr)e""'^ form the standard ba- 
sis for billiard problems, k is the momentum, E — k^, 
and Cm = Cl„i ^'"^ independent random variables with 
Gaussian distribution. Only numerical calculations of 
this distribution have been performed in Q| . 

The purpose of this letter is to demonstrate that nodal 
domains of random functions (and, consequently, 
wave functions of generic chaotic systems @1 , ) can be 
described by a simple percolation-like model where all 
interesting quantities can be calculated analytically. The 
model permits also to apply ideas and methods devel- 
oped within the percolation theory to the field of quan- 
tum chaos. 

To understand how the nodal domains look like we give 
in Fig. ^ their picture for random function (^ with k = 
100. The figure corresponds to a square window of the 
size L = 4 which contains 907 connected nodal domains. 
The largest of them and the largest of the domains which 
do not touch the boundary are highlighted. 



Our first step is to calculate the mean number of zeros 
of random functions (|l|) along a straight line (say the 
vertical one). This can be achieved by noting that, if the 
size in y-direction is Ly, the approximate quantization 
condition reads kyLy « nm where m is an integer and 
ky is the mean square momentum along the y-axis, k^ — 
k^ /2. Therefore when x is fixed 
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FIG. 1. Nodal domains for a realization of random func- 
tion (|^) with k = 100. Two different connected domains are 
highlighted. 

The rigorous derivation of this relation using the 
method of |^ will be given elsewhere |^. The same an- 
swer, of course, can be obtained for the mean density of 
nodal lines along any other straight lines ^ . 

These simple arguments show that the mean density of 
intersections of any straight line with nodal lines of ran- 
dom functions (|l|) is the same as in ^ and nodal lines of 
these random functions can in the mean be considered as 
forming an approximate rectangular lattice whose total 
number of sites is asymptotically 

k^A 
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-N{E), 
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where A is the area of the billiard and N — AE/Air is 
the mean number of levels (for billiards) with energy less 
than E. 
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But this simple picture can be valid only in the mean. 
When ^'(x, y) = ^{x, y)+5^{x, y) where function ^{x, y) 
has a crossing of nodal lines as in Fig. ||a, the addition 
of a small correction (5^'(x, y) will change, in general, the 
true crossing to one of two possible avoided crossings as 
in Figs. ||b and 
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FIG. 2. a) True nodal crossing, b) and c) Avoided nodal 
crossings. 

Consequently, one can conjecture that the distribution 
of nodal domains for random functions is the same as 
for the following random percolation-like process. Let us 
consider a rectangular lattice with the total number of 
sites Ntot = 2N{E)/tt as in (|). Each line crossing with 
probability 1/2 is changed either to the avoided crossing 
as m Fig. |b or to the one as m Fig. IId. These rules 
give a well defined random percolation-like process. One 
realization of such a process is presented in Fig. |^. 
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FIG. 3. A realization of random percolation-like process. 
Plus and minus form two dual lattices. Solid and dashed 
lines indicate graphs for respectively negative and positive 
dual lattices. 

The original lattice gives rise to two dual lattices called 
below positive and negative whose vertices are in the cen- 
ters of regions where our function is positive or negative 
(see Fig. ^ and whose size, a, coincides with de Broglie 
wave length (cf. (||)): a ~ 27r/k. Any realization of the 
above mentioned random process uniquely defines two 
graphs on these lattices, (which we call also positive and 
negative) with the following properties (i) their vertices 
coincide with the vertices of the corresponding lattice, 
(ii) their edges join together the connected components 
of this lattice. (A point is also a component of the graph.) 



One can choose arbitrarily a graph on one lattice (say 
negative) and any of such graphs will correspond to an 
allowed realization and vice versa. Therefore our ran- 
dom process is determines mostly by the bond percola- 
tion model on one of dual lattices (see e.g. [^) where with 
probability 1/2 one connects 2 nearby sites by a bond. 

The number of connected nodal domains coincides with 
the sum of the numbers of different components of both 
positive and negative graphs. As in [Q we first are inter- 
ested in the distribution of these numbers. To compute 
this quantity (unusual for the percolation) it is conve- 
nient to connect this model with the Potts model (see 
e.g. 0]) similarly as it was done in j|] for a slightly dif- 
ferent problem. 

Let n± be the numbers of connected components of 
positive and negative graphs. The generating function of 
their sum is 
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realizations 



(4) 



where variable x plays the role of the fugacity. 

The negative and positive graphs, by construction, are 
dual to each other Q and their properties are interre- 
lated. In particular (see e.g. p. 242) n+ C_ + 1 
where C_ is the number of independent circuits on the 
negative (dual to the positive) graph. According to the 
Euler relation this quantity can be expressed as follows 
C_ — b- + n_ — N- where b- is the number of bonds, 
ri_ is the number of connected components and is 
the number of vertices of the negative graph. 

These relations permit to express the generation func- 
tion (^) through the properties of only negative graph, 
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where we take into account that N- equals the total num- 
ber of sites of negative lattice, = Ntot/2. 

But this quantity is directly connected with the parti- 
tion sum of the Potts model ||^, The later can be 
defined for an arbitrary graph by the formal sum 



Z 



Pott, 
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where the summation is performed over all graphs, G, 
which cover the original graph. 6(G) is the number of 
bonds of this graph, n(G) is its number of connected com- 
ponents, q is the number of states of the Potts model, 
V — — 1 is a parameter related with the inverse tem- 
perature K . 

Comparing (H) and (^ one gets 



Z{x) 



■iPotts 



{x,x ). 



(7) 



The last sum corresponds to the Potts model in the crit- 
ical point = q; for large rectangular lattice and q < A 
it was computed analytically [fol 
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lim — In Zpottsix,x'^) =logx + f{x), (8) 



where 



dt , sinh(7r-u)t 
f(x) = / — tanh^t- ^ ^' 



t 



sinh nt 



(9) 



and the parameter /i (0 < /i < 7r/2) is related to the 
fugacity x as follows: cos^ — x/2. 

The expansion of this integral into the series on z gives 
the number of nodal domains with a fixed number of 
components, iV„. With exponential accuracy 



^7V„z" = zexp(7V,/(z)). 



(10) 



agree with the predictions ( p^ and (JIJ). Small errors 
may also be related to finite size effect as the infinite 
sum in (|l|) has been cut at |m| < 2kL where L = 4 is the 
side of the chosen square window. 

The close relation of our model with the bond per- 
colation model at the critical point p = \/2 permits to 
apply all results of percolation theory (see e.g. Q and 
references therein) for description of nodal domains of 
random functions. 
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We are interested in the behavior of iV„ for large n near 
the maximum of iV„ with fixed number of sites. One has 
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where the integration is performed over a contour around 
zero. Assuming that n is large and using the saddle point 
method one obtains Nn oc exp <i>(n, Zc), where $(n, z) = 
Nsf{z) — nlnz and the saddle point Zc is determined 
from the equation d^{n, z) / dz\z=zc = 0. The maximum 
of Nn corresponds to n — fi for which $(n, Zc) is maximal. 
Simple calculation shows that it appears when Zc = 1 and 
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Expanding $(n, z^) near n = fi up to the second order 
and computing all necessary integrals (the details will 
be given elsewhere [Q) one finds that the total number 
nodal domains in the lattice with Ntot = 2N{E)/-k sites 
(cf. (|^)) near the maximum has Gaussian distribution 
where the mean number of nodal domains, n, and their 
variance, cr^, are related to the mean number of levels in 
the following way 



n = N{E) 
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These relations are the main analytical result of this note. 
They demonstrate that the distribution of nodal domains 
for random functions is universal (Gaussian) and depends 
only on the mean number of levels exactly as it is con- 
jectured in ||l| with explicitly calculable constants. 

In Fig. ^ we present the comparison between the nu- 
merical calculation of the mean value and the variance for 
the random functions (|l]) for different values of k plotted 
as functions of = N{E). The dots and squares repre- 
sent the ratios n{N)/N and a'^{N)/N averaged over 100 
realizations. Within statistical errors numerical results 
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FIG. 4. Mean values of the nodal domains (dots) and their 
variances (squares) for random functions divided by N versus 
TV. The solid and dashed horizontal lines represent theoretical 
asymptotic values (^3|) and ( |l^ respectively. 

In particular the percolation theory predicts that the 
distribution of the areas, s, of clusters (= connected 
nodal domains), n(s), should have power behavior 



n(s) oc s 



(15) 



where the Fisher exponent r — 187/91 |^ p. 52. 

In Fig ^ we present the results of numerical calculations 
of this quantity for random functions (|l|) with A; = 115 
which are in a rather good agreement with the perco- 
lation theory prediction. In this figure the j/-axis rep- 
resents the number of nodal domains divided by N(E) 
and the area along the x-axis is measured in the unit of 
Smin = ■"■(ji/fc)^ where ji is the first zero of the Bessel 
function, Jo(ji) — 0, which according to the Rayleigh in- 
equality |]ll| is the smallest possible (with fixed k) area. 
After such scaling the results for random functions with 
different k are practically superimposed. The existence of 
discrete set of smallest possible areas leads to pronounced 
oscillations at small s in Fig |[ 

Another interesting quantity is the fractal dimension of 
the nodal domains. In our percolation model it coincides 
with the fractal dimension of critical percolation clusters 
which is known to be equal to I? = 91/48 B p. 52. 
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FIG. 5. Distribution of nodal domain areas. The solid line 
has the slope r = 187/91 predicted by the percolation theory. 

To find numerically fractal dimension of a domain it is 
convenient to put it on a grid of squares of side R and 
count the number of crossing of the region with the grid. 
When a <C i? <C / where / is the size of the domain and 
a is the size of the mean lattice: a — 27r/fc, one expects 



n oc R 



-D 



(16) 



and the exponent D is the fractal dimension. 

In Fig. ^ we present numerical verification of this re- 
lation for the two nodal domains with k — 100 high- 
lighted in Fig. [| and for the largest cluster in the pro- 
posed percolation-like model with number of sites given 
by (III). It is clearly seen that both domains have close 
fractal dimensions which agree well with simulations in 
our percolation-like model and the percolation theory 
prediction. 

To summarize, we developed a simple percolation-like 
model to describe the nodal domains for random func- 
tions. Its main advantage is that all relevant quantities 
can be computed analytically. By using the relations 
with the Potts model we demonstrated that nodal do- 
mains are distributed according to Gaussian distribution 
whose mean value and variance are proportional to the 
mean staircase function with explicitly calculated param- 
eters (|3|) and ( |l4|) . Our results clearly indicate that the 
distribution of nodal domains for random functions is in 
the same universality class as critical bond percolation 
which permits us to predict different critical exponents 
like the Fisher exponent for the distribution of the nodal 
domain areas ( |l5|) and its fractal dimension (|6|). 

Many different generalization of the model considered 
are possible. We mention only the possibility to use non- 
critical percolation model for description of level domains 
of random functions, ^'(x, y) — e, with e 7^ 0. 
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FIG. 6. Number of intersections of nodal domains with 
square grid of size R. Dotted vertical line indicates the mean 
lattice size. Circles and squares correspond respectively to 
large and small highlighted domain in Fig. ^. The dashed 
line: results of numerical calculations for the largest cluster 
in the percolation-like model. The solid line: percolation the- 
ory prediction with the exponent D — 91/48. 
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